Three body problem in a dilute Bose-Einstein condensate 
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We derive the explicit three body contact potential for a dilute condensed Bose gas from micro- 
scopic theory. The three body coupling constant exhibits the general form predicted by T. T. Wu 
[Phys. Rev. 113, 1390 (1959)] and is determined in terms of the amplitudes of two and three body 
collisions in vacuum. In the present form the coupling constant becomes accessible to quantita- 
tive studies which should provide the crucial link between few body collisions and the stability of 
condensates with attractive two body forces. 
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Despite their diluteness atomic Bose Einstein conden- 
sates (BECs) are significantly influenced by inter-atomic 
forces. The binary potential enters their properties 
mainly through two body low energy scattering observ- 
ables described by the s wave scattering length a. Recent 
experiments have provided the opportunity to tune 
the scattering length from the repulsive (a > 0) through 
the ideal gas (a ^ 0) to the attractive (a < 0) regime 
where the gas finally becomes unstable with respect to 
collapse |§, |j. With this new technique one can thus 
access a wide range of dilute gas parameters and study 
phenomena related to few body collisions in the gas. 

Virtually all properties of BECs at zero tempera- 
ture are described by a nonlinear Schrodinger equation 
(NLSE) i 



■U^ T 



fflB* + .g2|*P* + 331*1^* 
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which, in this extended form, depends on the Hamilto- 
nian of a single trapped atom as well as the two and 
three body coupling constants §2 and g^, respectively. 
Three body scattering obviously becomes significant at 
large scattering lengths, through the increase of the dilute 
gas parameter ij — y/\W^]a\^, but also in the ideal gas 
regime. Small scattering lengths of only a few A, e.g., are 
needed to produce stable condensates with attractive in- 
teractions and several thousands of atoms in present day 
atom traps. Two body scattering alone, however, seems 
not sufficient to explain the experimental criterion for 
stability against collapse of *^Rb |^ and a fit of Eq. (|l|) to 
the data gives an attractive three body interaction with 
very reasonable coupling constants jRegal/ft of the order 
of 10~^^ cm^/s 1^. In this way the corresponding three 
body mean field shift of the energy density, Regsj'I'l^/S, 
should become directly accessible to experiments. 

The determination of the three body coupling constant 
in a dilute BEC, and its associated mean field energy, 
have a long history of theoretical research in many body 
physics Already in 1959 Wu ||^ discovered the gen- 
eral form 33 — 16nh'^a^{4:7r — 3\/3) ln{Cri'^)/m for a Bose 
gas of hard spheres. The constant C in the argument of 
the logarithm, however, remained undetermined and only 



recently was significant progress made. In 1999 Braaten 
and Nieto used effective field theory to determine 
Wu's general result by taking into account low energy 
three body scattering. They estimated the argument of 
the logarithm on the basis of general assumptions on the 
length scales set by the inter-atomic potential. Braaten 
et al. 1^ also studied the role of parameters of the po- 
tential apart from the scattering length, through fits to 
many body Monte Carlo simulations [|l^, as well as the 
properties of the effective three body coupling constant 
when the scattering length is large in comparison to the 
range of the potential (ll| . 



To the best of our knowledge, the three body coupling 
constant has not been derived from a microscopic de- 
scription of the gas and expressed in terms of the binary 
potential. In this letter we provide a microscopic theory 
of three body collisions in a BEC and determine g2 and 
g3 to second order in the dilute gas parameter rj. The 
derivation is based on the exact quantum evolution of 
correlation functions, i.e. quantum expectation values of 
products of single mode annihilation and creation oper- 
ators of the form (ap), (ap^Up^), (aqOp) etc.. References 
]l2| , provide a general approach to transform the ex- 
act infinite hierarchy of dynamic equations for correlation 
functions into a more favorable form that allows for a sys- 
tematic truncation in accordance with Wick's theorem. 
The resulting set of dynamic equations for non commu- 
tative cumulants lends itself to iteration. We shall 
apply this cumulant approach here to express the exact 
dynamic equation for the condensate wave function of a 
uniform gas through nonlinear terms in ^ and ^* and 
determine g2 and g^ in equilibrium. 



According to Ref. |l3| the exact quantum evolution of 
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order cumulant A which is associated with the correlation 
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formal solutions then read: 



FIG. 1: Jacobi coordinates of three atoms. The spatial coordi- 
nates are chosen as the vector from atom 1 to atom 2 (ri2) and 
the vector from the center of mass of atoms 1 and 2 to atom 3 
(p). The conjugate Jacobi momenta are p — (p2 — pi)/2 and 
q = (2p3 — pi — P2)/3, respectively, where Pi,P2,P3 denote 
the momenta of the single atoms. 



gas limit. The dynamic equation for \1/ then reads ^ 
1 




dqdpT/(pi)A(q,p,t) 



dpV{p) $(p,t) + V2^ 5(p)vl/2(t) M,*{t), (2) 



where the integration variables are chosen as center of 
mass and relative atomic momenta q = pi + P2 and p = 
(P2 — Pi)/2, respectively, and V{p) = J drl/(r) exp(— ip • 
r/h)/{\/2nh)^ is the Fourier transform of the potential 
into momentum space. Additional terms [p^ that are 
proportional to the density matrix of the non condensed 



fraction, Fpq — 



— (aq)(ap), neither contribute to 
f/2 nor to (73 and have been omitted in Eq. (||). 

The explicit dependence of Eq. (||) on higher order cu- 
mulants can be eliminated successively by solving their 
dynamic equations formally and inserting the solutions 
into Eq. (||) iQ. With this iterative procedure Eq. (||) 
can be expressed exactly through nonlinear terms in ^P, 
"if* and r up to a given order and a remainder which 
contains only higher order cumulants. To obtain the cou- 
pling constants 52 and 53 we shall determine all the non- 
linear terms in ^E* and 4"* up to the fifth order. 

The formal solutions of the dynamic equations of $ 
and A p3| involve the triple function x, i.e. the third 
order cumulant associated with the correlation func- 
tion (flpjap^apj), and the Green's function correspond- 
ing to two atoms described by the two body Hamil- 
tonian H2B = —^^r + ^(r) in their center of mass 
frame. The relevant retarded Green's function is given 
by |l| G{t) = :^9{t)cjip{-iH2Bt/h). We shall use 
G{t, p,p') = {p\G{t)S\p') in momentum space, where 
S projects onto the symmetrized relative momenta that 
correspond to two identical Bose atoms. At the initial 
time ^0 the gas is assumed to be prepared in a (ideal 
gas) coherent state of the lowest energy mode IT^. The 



^P,t)= J^^dT I dp'G{t-T,p,p')V{p')^'{T) 

+ dr j dc{dp'G{t-T,p,p'2)[V{p',) + V{c{)] 



X(q',p',r) + 2V2^ 5(q')*(p',T)vI/(r) ^*{t), (3) 



and 



A(q, P. = ^ j dq'dp'Git - r, p, p^ - q/2) 
X X(q', p', r) + 2V2^'<5(q')<i>(p', t)*(t) 
X [Vip[ -q) + y(q' - q)]$*(q, T)e^^^(*-^)/'* 

dr j dq'dp'G(t-r,p,p')x(-q,p',T)T^(q + q') 



$(q',r) + v^2^ (5(q')^'^(T; 



(4) 



where p'l = q'/2 — p' and p'2 = q'/2 + p'- The form of 
Eqs. (H) and (^) is exact, save that terms that neither 
contribute to 32 nor to (73 are not given. The arguments 
of the triple function x for a uniform gas are chosen here 
as the three body Jacobi momenta in Fig. ^ 

The dynamic equation for x involves the Hamiltonian 
of three atoms interacting pairwise. 



+ Ax,,) 



(5) 



where V^b = V{^2z) + l^(ri3) + V{v^2) (see Fig. |l|). The 
corresponding symmetrized retarded Green's function in 
the three body center of mass frame is denoted here as 
G3B(i, q, p, q', p') in the representation through the Ja- 
cobi momenta in Fig. |l]. The only relevant contribution 
of the triple function reads 

x(q,P,<)=6^ dT J dq'dp'G'3B(t - r,q,p,q',p^) 

xF(q')$(p',r)vI/(r), (6) 

where P2 — q'/2 + p'. 

All nonlinear contributions in '5 and up to the fifth 
order are obtained by inserting Eqs. (||) and (^) into 
the right hand side of Eq. (|^) and then replacing x by 
Eq. (|^) and $ by its leading first contribution on the 
right hand side of Eq. (0). Equation (^ then becomes a 
non-Markovian NLSE [Oj with the leading (two body) 
contribution 



d f°° 
ih—-^{t)^{2Tihf / dr T(t-T,0,0)*2(r)^'*(t) 
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where T{t) = V6{t) + VG{t)V is the Fourier transform 
of the energy dependent two body T matrix [nsl. This 
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time dependent transition matrix is sharply peaked at 
t = with a width determined by the two body coUi- 
sional duration. The fifth order nonhnear terms in 5* 
and (three body contributions), indicated by the dots 
on the right hand side of Eq. (|^) , involve up to four time 
integrals which have a similar form to the leading term. 

To determine the equilibrium properties of the gas we 
shall assume first that the scattering length is positive 
and the potential does not support any bound clusters of 
atoms. In the zero temperature equilibrium the ansatz 
\Sf{t) = ^/ exp{-iij.t/h) then transforms the NLS Eq. (|^), 
including the three body contributions, into its station- 
ary form in the limit t ~ to oo |17| : 



(8) 



The two body contribution to Eq. (|^) is directly ob- 
tained from Eq. (Q) through the Fourier transform of 
T{t—T, 0, 0) at the energy 2/i which yields the two body T 
matrix T(2//, 0, 0) (18). The low energy expansion H, |||] 
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shows that the leading contribution to /i is the mean field 
chemical potential /ig = 47r?i^a|^p/m. The three body 
contributions to the right hand side of Eq. (^) (not shown 
explicitly) can also be arranged in such a way that all 
bare two body potentials are connected to two and three 
body T matrices. The general form of the real part of 
Eq. (^ is then obtained as (77 = ^j^'pa^): 
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where corrections smaller than ^iqT]^ have been neglected. 
The first term on the right hand side of Eq. ( [lo| ) is the 
two body part. All three body contributions to Eq. (H) 
are at least of third order in the two body T matrix. The 
real parts of the third and fourth order terms exhibit a 
logarithmic divergence at low energies /i — > as indicated 
in the second term on the right hand side of Eq. (|l0|). The 
respective explicit expressions read 



In 



P 



(3) 



-Re 



X ^[r(2M,p,0)] 

a=±l 



rfpr(2A.,p,o) 

2/i — /m + iO 



4m ' " 2 ' 2> 



[T(2a*,p,0)] 
2T(2Ai,0,p) 



2/i — /m — iQ 

, T(/.+ i^,aH £) 
2/i — jm 



4in ' " 2 ' 2) 



2fj, — p"^ /m + iO 



(11) 



and 
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Q,/3 = ±l 



T(- 
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The logarithmic behavior of the real parts of the integrals 
in Eqs. (pi] ) and (|l2|) was determined with the methods 
in Ref. | pO| . The momentum units p'^' and p^''^ are found 
using the regular parts of the integrals in the limit /i — *■ 0. 

The last contribution to the right hand side of Eq. ( |l0| ) 
involves the T matrix for three asymptotically free incom- 
ing and outgoing atoms, T^b{E) = V3B + X^3BG3B(£^)^3B■ 
(5) 

The regular part of T3B , denoted here by Tgg , consists of 
all multiple scattering contributions to Tsb from the fifth 
order upward [|l|. T^^^{E) can be evaluated at ii^ = 
po| and is the only contribution to 173 that depends on 
the complete zero energy three body scattering state. 
Equation (ph also exhibits an imaginary part 



Im/i = HQ 



m2/i 



y/m2iia 



Iff- +0{rf Inry) 



(13) 



whose first leading term stems from Eq. (g). The second 
term is the leading contribution to the imaginary part 
of the integral in Eq. (11). The imaginary part of /i 
describes the exchange of atoms between the condensed 
and thermal fraction jist and is thus required to vanish 
exactly in equilibrium as soon as all few body scattering 
contributions are taken into account. The right hand 
sides of Eqs. (^0|) and (^3|) change by an amount smaller 
than iiQif if ^ is replaced by its leading contribution fiQ. 
The square root terms in Eq. (13[) then, indeed, cancel 
each other and Eqs. 




.92 =(27rft)3ReT(2/io,0,0), 
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A trapped dilute BEC can usually be divided into re- 
gions with an extent much larger than the range of the 
inter-atomic forces but still much smaller than the scale 
of the spatial variation of the condensate wave function 
\l/(x). In these regions the gas can be considered as lo- 
cally homogeneous and the coupling constants in Eq. ([l^ ) 
remain valid when /io is replaced by the local chemical 
potential fJ-o{^) — 47rft^a|^'(x)p/TO. 

A realistic molecular potential V{r) usually supports 
two body bound states. The coupling constant 33 then 
acquires an imaginary part through t!^^ . According to 
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three body collision theory this imaginary part corre- 
sponds to the threshold-less formation of dimer molecules 
and was determined here via an optical theorem with 
the methods of Ref. ||2^. Wc thus obtained 2lT[ig^/h = 
—K3/6 where is the three body recombination rate 
constant for an ultra-cold non condensed Bose gas. This 
was determined in Ref. ||2^ using classical probability ar- 
guments, by assuming that three atoms are lost in each 
event. The factor of 1/6 expresses the fact that all three 
condensed atoms share the same quantum state |Q and 
was verified in subsequent experiments . In the pres- 
ence of bound states the condensate can only assume a 
metastable state. The change of its density, however, is 
negligible on time scales comparable to collisional dura- 
tions so that the present results remain valid. 

If a is negative but the trapped gas is still sta- 
ble with respect to collapse the above analysis slightly 
modifies: /io(x) then becomes negative and instead of 
Eq. ( p^ ) one obtains real contributions, /iq ^/rrl2\f^\a / h 
and ^o87rft|5'pa^/-\/TO2|^|, through Eq. (||) and the in- 
tegral in Eq. (pi]), respectively. The contribution of 
Eq. (pi]), with \^\ replaced by |/zo|, revises the two 
body couphng constant to 32 = (27rfi,)'^r(2/io, 0, 0) — 
ATTha^y/2]LLo\/rn. This leads to 52 = ^fa[l + 0{ri^)] 
by Eq. (||). If the gas collapses |^ the equilibrium as- 
sumption as well as the subsequent concept of two and 
three body coupling constants break down. 

When a three body bound state emerges close to the 
zero energy threshold (73 is dominated by the residue of 
the corresponding pole through T^^ The real part 

of 53 can then assume all positive and negative values 
in complete analogy to the two body case. This phe- 
nomenon was discussed in Ref. [Tl|] in connection with 
the Efimov effect in the three body energy spectrum [ p6[ . 

The three body coupling constant in Eq. (|l4) has the 
form predicted by Wu 0] and Braaten et al. ||, ^ and 
can be determined quantitatively by solving the two and 
three body scattering problem in vacuum [ p7| . The ex- 
plicit representation in Eq. ( p^ also allows to directly 
relate all our current knowledge about two and three 
body transition amplitudes to properties of the gas. The 
present theory is valid for all (non singular) molecular 
potentials and avoids divergences. Convergence at high 
energies has been achieved by taking into account the ex- 
act binary potential which enters the coupling constants 
through two and three body T matrices. At low energies 
the arguments of the infrared divergent contributions to 
53 have been determined to multiples of fiQ which is the 
energy scale that corresponds to the healing length 

In summary, we have derived the explicit three body 
contact potential for a dilute condensed Bose gas from a 
microscopic many body theory. The respective coupling 
constant, g^, that reflects general equilibrium properties 
of an interacting BEC has been expressed in terms of the 
binary potential through two and three body transition 



matrices. Its explicit form allows for a determination of 
the magnitude of from few body physics. We have 
derived the heuristic three body recombination rate con- 
stant of Ref. as well as phenomena related to resonant 
three body bound states. 

I am greatly indebted to Thomas Gasenzer, Simon 
Gardiner, Robert Roth, and Keith Burnett for valuable 
discussions. This research was supported by the Alexan- 
der von Humboldt Foundation and the UK EPSRC. 
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